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We derive the large deviation function, which provides the large-time full counting 
statistics for the charge transfer, in the non-equilibrium steady state of the resonant-level 
model. The general form of this function in free fermion models, in terms of transmission 
coefficients, was proposed by Levitov and Lesovik in 1993 using a particular measurement 
set-up involving an interacting spin. It was later suggested to hold as well for a proper 
quantum mechanical measurement of the transferred charge. We give a precise proof of 
both statements in the resonant-level model. We first give a full description of the model 
and its steady state. That is, we explain how the decoupled system prepared with a charge 
differential evolves, with the impurity coupling, towards the Hershfield non-equilibrium 
density matrix, in the sense of averages of finitely-supported operators. We describe how 
this holds both for the usual resonant-level model with a point-like impurity, and for a 
regularised model with an impurity spread on a finite region, shedding light on subtleties 
associated to the point-like impurity. We then prove Levitov-Lesovik formula by recasting 
the problem into calculating averages of finitely-supported operators. 
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1 Introduction 



Full counting statistics in electronic transport refers to the probability distribution func- 
tion (p.d.f.) of charges transferred through quantum wires. In the simplest setting, one 
may have in mind (or for real) a quantum dot coupled to (at least) two electronic reser- 
voirs which are maintained at different potentials (and different temperatures) so that an 
electronic current flows through the dot from one reservoir to the other. The system is 
by construction out of equilibrium. Since quantum effects dominate, one has to describe 
the reservoirs - and, of course, the dot - quantum mechanically. 

Let q be the charge transferred during time t. Its p.d.f. is fully characterised by the 
generating function P{X,t), 
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with Ilq the probability of observing a transferred charge q. A formula for P(A, t) at large 
time has been given by Levitov and Lesovik |4]. However, one has to specify how the 
transferred charges are actually measured. 

A first way to define the transferred charge during a time duration t consists in first 
measuring at time the charge, say, in one of the reservoirs, letting the system evolve, 
and then measuring again the charge at a later time t. Let Q be the charge operator 
and Pq the projector onto the Q-eigenspace with eigenvalue q. If po is the normalised 
system density matrix at time 0, basic principles of quantum mechanics tell us that the 
generating function for the transferred charges is: 

P(A,t) = 5]e^^''Tr(p,„+,t/o,P,„ponot^o^A+9) (1-1) 

where Uo-^t is the system evolution operator from time to time t. The sum over go 
corresponds to the sum over the values of the charge measured at time 0, while that over 
q + qo corresponds to the charge measured at time t. Note that this formula involves two 
projectors Pg^ and Pg+go- 

Definition (11.11) is actually not the way Levitov-Lesovik [1] chose to specify the charge 
transferred. Instead they define the measured transferred charge by modeling a measure- 
ment device made of a quantum spin coupled to the electronic current flowing through 
the quantum wire. A precise description of the measurement system is not necessary for 
our discussion but may be found in [3] . We only recall that it leads to another definition 
for the generating function: 

P(A,t) = 5^e'^''Tr(Pg,+gf/o,Pg„pof/oV) (1-2) 

9,90 

with, again, po the system density matrix and Uq-i is the evolution operator. This formula 
also involves the two projectors Pg^ and Pg+go but it differs form the previous ones by 
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the fact that there is only one insertion of the projector Pg^ and not two as in Eq. fll.ll) . 
The two formulae would coincide if the density matrix happens to commute with the 
charge operator, and thus with the projectors Pg. But this is generically not the case, and 
certainly not if the density matrix represents a steady state where the charge Q is flowing. 
None of these formula correspond to log Tr(e*'*'^*^*^'^'' p) which would be the analogue of 
the classical formula. 

The universal character of the p.d.f. encoded in fll.ll) only emerges at large time 
duration t and in a stationary regime, that is only when the system "reservoirs plus 
dot" is in a stationary state described by a (non-equilibrium) density matrix Pstat- Let 
-Pstat(A,i) be the corresponding generating function. We are interested in evaluating the 
large deviation function F{X) defined by 

F(A) = limrMogPstat(A,t). (1.3) 



The Levitov-Lesovik formula [1] is an expression for -F(A) associated to the generating 
function ( II. 2p . in the case where the electrons propagate in the quantum dot through one 
channel only. It reads: 



F(A)^/^.og(l.T( 



(to) n^{co){n2{co)-l){l-e'^)+n2{co)Mco)-l){l-e-'^) ) (1.4) 



with T{uj) the transmission coefficient at energy oo and 



nj{u) 
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(1.5) 



the Fermi occupation numbers of the two leads. When the two reservoirs are at identical 
temperature, Ti = T2 = T, the Levitov-Lesovik large deviation function satisfies the 
fiuctuation relation [121 M- = F{—\ — iV/T) with V = Hi — ^2 the potential 

difference. 

Of course, there is already quite a number of papers dealing with the Levitov-Lesovik 
formula. The original paper [1] as well as the very nice review by Levitov [5] deal with 
the definition (11. 2p of the transferred charge motivated by a model for the quantum 
measure apparatus. Recall that definition (II. 2p does not coincide with the von Neumann 
prescription for the measure of the quantum charge at different time intervals. Levitov- 
Lesovik [U inj derivation is based on a scattering approach. References [6], [7] also deal with 
definition (11.21) of the transferred charges either using a microscopic formulation [7] or a 
scattering approach combined with free fermion technology [B]. References [S], [H] start 
from the von Neumann definition of the measured transferred charge but Ref. [8] assumes 
that the charge operator commutes with the density matrix and Ref. [9] uses a quantum 
master equation approximation. Up to our knowledge, the approach to the stationary 
regime in relation with the full counting statistics does not seem to have been analysed 
in the previous literature, except partially in Ref. [TU]. 
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The aim of this paper is to provide a complete derivation of the Levitov-Lesovik 
formula in the framework of the resonant-level model (whose definition is given below), 
one of the simplest models for a quantum dot. Although the formula is already known 
for some time, the derivation in the framework of the resonant-level model contains a few 
subtleties which have never been thoroughly sorted out. We introduce a regularisation 
of the model in which the impurity is spread over a small distance. First we prove that 
the system reaches a stationary state, which requires taking large size and large initial 
preparation time limits in the appropriate order. This limit actually exists only for local 
operators (as was suggested in the context of the Kondo model [IQ]), that is, operators 
localised around the quantum dot. The charge operator is not a local operator, since it 
measures the charge in a reservoir, so that it does not admit large size stationary limit. We 
however show that the operator involved in the transferred charge p.d.f. is nevertheless a 
local operator and that it admits a stationary limit. We then use free fermion techniques 
to evaluate the large deviation function. We also show that the two ways of defining the 
measured transferred charge, fll.ip and fll.2p . lead at large time and in the stationary 
regime to identical large deviation functions. 

Our approach is strongly based on local-operator methods, and in that differs signifi- 
cantly from previous approaches to derive the Levitov-Lesovik formula. The idea of local 
operators is crucial both for the explicit calculations, and for showing the existence of 
the stationary limit. On the other hand, techniques using scattering states may be more 
amenable to interacting, integrable models, for instance. We did not make any attempt 
to connect our derivation to the scattering-state approach, and hope to come back to this 
in the future. 

The paper is organised as follows. In Section O we define the resonant-level model 
taking care of the subtleties introduced by an a priori point-like impurity, we discuss how 
to describe the stationary state using density matrices of the Hershfield type, and we show 
that the stationary state is reached in the appropriate large-time limit. In Section |3l we 
prove the Levitov-Lesovik formula for the resonant-level model, using local-operator and 
free-fermion techniques. 

Acknowledgment 
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2 Non-equilibrium stationary state in the resonant- 
level model 

2.1 General discussion: model and the stationary state 

The resonant- level model describes a quantum dot coupled to two electronic reservoirs. 
Each of the reservoirs is modeled by massless fermions propagating on a line segment 
[0,L]; after unfolding, this becomes chiral fermions on the interval [—L,L] with periodic 
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boundary conditions. The dot is modeled as simply as possible by a unique level which 
can be either empty or occupied. The coupling is such that electrons may hop between the 
dot and the reservoirs with identical amplitudes for hoping in either of the two reservoirs. 
The Hamiltonian we consider is H = Hi + H2 + Hd + i^int with Hj the free fermion 
Hamiltonians describing the reservoir dynamics, Hd the dot Hamiltonian, and //mt the 
dot-reservoir coupling. Explicitly [Ij: 

H = -iVF J ^ dxi4d^iJi + + ^((^1(0) + ^2(0))rf + d\MO) + Mm + e dU. 

(2.6) 

where is the fermion creation operator on the dot, {S, d} = 1, e the dot energy, ipj{x), 
j = 1,2, are the fermionic operators of the reservoirs, {■?/'] (x), ipk{y)} = Sj;kS{x — y), with 
periodic boundary condition ipj{L) = tpj{—L), vp is the Fermi velocity and r the hopping 
amplitude. 

The system is placed out-of-equilibrium by preparing, say at early time to < 0, the two 
reservoirs at different temperatures, Ti and T2, and different chemical potentials /ii and ^2- 
The initial system density matrix is pd ® pth with pth the reservoir thermal density matrix, 
with temperature Tj and chemical potential pj, and pd the dot density matrix. Explicitly, 
p^^ = pi (g) p2 with pj = Z-^e-l^^^^^-''^^^\ Tr(pj) = 1, where Nj = J dxip]{x)ipj{x) 
and Hj = —ivp J_j^dx i/jjdxi/jj (here Pj = 1/Tj). Upon shifting energies, we may take 
PiPi + /32p2 = without loss of generality; we will assume that this relation holds in the 
following. 

Charge and heat currents are going to flow between the reservoirs through the dot. 
They roughly propagate at velocity vp, so that a stationary regime is reached on a domain 
of scale £stat (this is a quantum correlation length, here proportional to r^) around the 
dot for large enough early time |to| and large L with 

L > Vp\to\ > istat- 

The reservoir size has to be large enough for the charge and heat flows not to bounce back 
on the reservoir boundaries (see Fig. [1]). 

Physics in the stationary domain is described by a stationary density matrix pstat of 
the coupled system which, according to the previous discussion, is given by infinite L and 
to limits in the following order: 

Pstat = hm lim Uo-to ( Pd ® Pth ) UL (2.7) 

to— >— 00 L— >oo 

with Ut-to = exp(— z(t — tf))H) the time-evolution operator from to to t. This limit exists 
(in the weak sense) as long as the density matrix is evaluated against operators localised 
in the stationary domain |x| < Igtat ^ ^fI^oI- 

The stationary state density matrix can be described using Hershfield's density matrix 
pT| (the following is a slight generalisation, with /3i 7^ in general): 

Pstat = e-^^-^^+^^ (2.8) 
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Figure 1: A (cartoon) picture of the two (unfolded) electron baths and the impurity, and 
of the evolution of the system. The figure on the left represents the initial density matrix: 
the upper electron bath is filled and the lower one is empty, due to the initial potential 
difference. The large arrow indicates time evolution. The figure on the right represents 
the density matrix after a time |to|, evolved without a potential difference, along with 
the instantaneous motion of the electrons. They fall from the filled electron bath to the 
empty one in order to establish their zero-potential configuration, and the region of fallen 
electrons expands with velocity vp- The steady-state limit is the one where the system is 
much larger than the region of fallen electrons, and where this region is much larger than 
the quantum correlation region around the impurity. 

where /3 = (/3i + /32)/2, V = — (32^2, W = f3i — (32, and the two operators Y and 
X are (generalisation of) Hershfield's operators describing the steady state. There are 
various ways one can formally define these operators. We will adopt that first used in [3] . 
It simply says that Y ~ (A'^i — N2)/2 and X ~ {Hi — H2)/2 for x < (in general, for x 
small enough), and that both Y and X are conserved by the dynamics. That is, we have 
[y, H] = [X, H] = and: 

[X,iIji{x)] = ^wpd^-^iix), forx<0; (2.9) 

[X,ip2{x)] = -livFdxip2{x), for x < 0. 

These specifications of the non-equilibrium density matrix apply to models defined as 
in Eq. (12. 6p with a localised impurity. We shall introduce a regularised version of the 
model in which the impurity is spread over a small distance. The above relations (12. 9p 
will then be valid only for x negative enough to be away from the impurity support, 
instead of all x < 0. Note that pstat describes a current flowing from lead 1 to lead 2. 
Physically, the conditions indicate that pstat is stationary and that pstat reproduces the 
initial density matrix at positions where the electrons have not yet interacted with the 
impurity, i.e. before they flow from one lead to another. We will see below that there 
is a unique solution (up to the addition of terms proportional to the identity operator) 
to these conditions. Although the density matrix has a form that resembles that of an 
equilibrium density matrix, as was argued in [3l [11] it is the fact that the operators Y 
and X are non-local which makes it a non-equilibrium density matrix. 

In the following, we set = 1. In Subsection 12. 2[ we recall and develop further the 
solution to the resonant-level model, both by using a "resolution" of the impurity, and 
by using an explicit spreading of the impurity, in order to take care of the ambiguities 



[r,V^i(x)] =--^i(x). 
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produced by the impurity-interaction terms in the Hamiltonian f l2.6p . In Subsection 12.3^ 
we prove that (12 .Qp has a unique solution, and we prove (12.70 . 



2.2 The operator solution to the resonant-level model 

Let us now recall the solution of the resonant-level model (the diagonalisation of the 
Hamiltonian). We use the even and odd combinations 

1 1 

in terms of which the Hamiltonian decouples. As usual, the approach is to use a parametri- 
sation of the operators ipe{x), ipo{x) and d in terms of mode operators, such that the 
canonical anti-commutation relations hold and such that the time evolution of the mode 
operators is simple. We will look at various ways of regularising the impurity term present 
in the Hamiltonian (hence we will have various time evolutions). 



2.2.1 Partial regularisation via a resolution of the impurity 

The solution to the odd part is trivial as it does not involve the impurity (see (12.170 below) , 
hence let us discuss the even part. Since the Hamiltonian is quadratic in fermion operators, 
the full solution can be straightforwardly obtained from the one-particle solution. The 
one-particle (even) Hilbert space is the direct sum L^(§)©C: a square-integrable function 
g{x) on the interval [—L,L] with g{—L) = g{L) representing the wave function of the 
even electron, and a complex number h representing the amplitude for the presence of the 
electron on the impurity. There are two equations resulting from the eigenvector equation 
for H with eigenvalue p: 

- id^g{x) + Th5{x) = pg (x) 

rg{0) + eh = ph. (2.10) 



Any solution to (I2.10p is such that the jump is g{0~^) — g{0~) = —irh. Yet, clearly, 
the second equation necessitates the knowledge of g{0), and this is also true for the 
evaluation of certain averages (e.g. the average of the momentum operator —idx)- A 
way of solving this problem is to use a resolution of the impurity point. We will replace 
ipe{0) in the Hamiltonian (12. 6p by the symmetric sum ipe{0) = ('^e(0'^) + '?/^e(0~))/2, with 
the understanding that 0"*^ are limits from the right/left towards 0. The delta-function 
potential becomes likewise a symmetric sum of delta functions, which, under integrations 
with functions in general discontinuous at a; = 0, takes their symmetric sums. Also, the 
second equation becomes well-defined, since g{0) is replaced by the symmetric sum of 
g{0^) and g{0~). This resolution is equivalent to defining g{0) as (^'(O"'') + g{0~))/2 for 
averages of operators with up to one derivative, and exactly reproduces the limit where 
a spreading of the impurity tends to the delta-function (see the discussion below of the 
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impurity spreading). We note that it is inconsistent to choose, as a resolution of ipei^), 
any non-symmetric sum of ilJe{0~^) and ipei^'), as noted in [2]; below we indeed show the 
universality of the chosen resolution: its independence from the way the impurity was 
spread. This resolution is not enough for an unambiguous definition of all averages, but 
it will be sufficient for our purposes. 

The Hilbert space is obtained as a completion of the space of solutions to the resolved 
equations f l2.10p . A solution is obtained by an appropriate Fourier-like transform. We 
show in Appendix O that the map (a^ : p) ^ g ®h given by 

(y'(x) = ep(x) Op, — ir /i = Op, (2-11) 

p p 

with 

eJx) := e^P" x | ^ for x < ,2.i2) 
^ ' t;p for X > ^ ^ 



and Vp and Wp defined by 



_ rV2 + .(p-e) _ 
"-=-rV2 + i(p-6)' ^- = ^'^-'=-rV2 + .(p-6)' ^^''^^ 

maps doubly infinite square- summable sequences (op : p) = {ap : e^*^^ = Vp) onto i^^(S) © 
C. It is a simple matter to see that sequences with only one non-zero element Cp yield 
solutions g = Cp, h = iwp/r to (12.101) . Note that the allowed values of p are obtained 
using the periodicity condition 

e^'P^ = Vp. (2.14) 

There is one eigenvalue p in each interval [ ^^^2L^^^ ; ''^^zl^''"^ ] except for the two intervals 
containing ±r^/2 which each contains two eigenvalues. We also show in Appendix O that 
the completion of the space of such solutions Cp © iwp/ r, under the usual inner product 

{gi® hi , g2® = dxgi{x)g2{x) + hih2, (2.15) 



is the Hilbert space L^(S) © C, and corresponds, under the map (12. lip , to the space of 
doubly infinite square-summable sequences (Cp : p). In particular, although every finite 
linear combination g ®h oi basis elements Cp © iWp/ r has the property that 

(7(0+)-(7(0-) = -zr/i, (2.16) 

this property does not survive completion: the function g and the number h may be 
chosen independently, in the sense that for any i^^(S) function g and h E C, there exists 
a square-summable sequence (op : p) that reproduces them under (12. lip . 
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From the one-particle problem we can solve for the fermionic operators dynamics. 
Since the set Cp © iwp/ r form a basis of eigenvectors of the one-particle Hamiltonian, this 
dynamics is diagonalised by introducing canonical operators and hp such that: 



—IT d(t) 



~2r 



-ipt 



2L 



V2n 

~2r 



e'P("-*) hp, 



(2.17) 



^Wpt 



-ipt 



with, again, quantification conditions e^*^'^ = 1 in the odd sector and e^*^^ = Vp in the 
even sector. From this, we see that the quantity 



Vp =: e 



i<j>P 



defined in (12.13P is the phase shift, and that the transmission coefficient is 

rV2 



Tip) = \wp/2\ = sin(0p/2) 



v/rV4+(p-6)2 



(2.18) 



(2.19) 



The inverse of the mode solution is 



(1 + ^)-^ 
^ 2Lr2^ 



dx 



IWr, 



-d 



2txt 



dx 



ipo{x)e 



~ipx 



(2.20) 
(2.21) 



The canonical commutation relations {'?/'|(x), ?/'j(?/)} = Si-j5{x — y) and {d\d} = 1 are 
equivalent to 

{bl,bp>} = ^6p,p>, {dl,dp>} = ^ (1 + J^y^^P'P'- 

The other anti-commutators vanish. 

The Hilbert space is the Fock space over these anti-commutation relations which makes 
the energy bounded from below. The vacuum |vac) is defined by 



a_p|vac 



) = ciplvac) = 6_p|vac) = 6],|vac) = W p <0. 



Its physical meaning is of course that the Fermi sea of negative-energy states has been 
filled. In the Hamiltonian, this infinite negative energy must be appropriately shifted 

away, so that 



H 



TT 

L 



J2p 



Wr. 



: I : ala^ 



2Lr2 J - P P 
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where the normal-ordering is with respect to the vacuum |vac). Note that thanks to 
fl2.16p . the operators tpe{x) and d are related to each other in any finite-energy state by 
the discontinuity relation 

iPe{0^) - M^") = -^rd (2.22) 

(this also holds in states containing infinitely many eigenstates of the Hamiltonian if 
the coefficients decrease exponentially with the energy) . The discontinuity relation comes 
from the fact that one set of modes is used to describe both ipe and d : this is a phenomenon 
usually referred to as hybridisation. 

We also need to take the large L limit. The density of eigenvalue of the one-particle 
Hamiltonian is uniform in this limit. Defining ap = limL_5.oo dp and bp = limL-s.00 bp, the 
Hamiltonian is simply H = J dpp^a^^ap + b^pbp), and the fermionic operators decompose 
as: 

i^eix^t) = / — =ep(x)e"*^*ap, 



/27r 

M^,t) = [ -^e^P(^-'Hp, (2.23) 
J y 271 

f diT} 

-irdit) = / —= Wpe'"^^^ ap. 



'2tx 

The canonical ant i- commutation relations give 

{alap>} = {blbp,} = 5{p-p'), (2.24) 
other anti-commutators vanishing. 



2.2.2 Pull regularisation via a spreading of the impurity 

An alternative way of dealing with the ambiguities of fl2.10p consists in regularising the 
Hamiltonian (12. 6p so that the impurity d interacts with the itinerant fermions on a neigh- 
bourhood of the origin. The regularised Hamiltonian is 

/h nL 
dx{i)%i)e + i^%i^o)+r / dx^^''\x){^l{x)d + Si)e{x))+ed^d, (2.25) 
-L J-L 

where if^'^^x) := a~^Lp{x/a) is a mollifier of the Dirac distribution. The function is real 
with compact support on M, normalised such that J^^dx (p{x) = 1. We consider a small 
enough so that the support of (^'"1 lies in the open interval {—L,L). We will denote by 
la = [la^-^a] closcd interval delimited by the maximum/minimum values of the 
support of cp^"''^. We shall be interested in the limit a — )■ 0"*" in which (p^°'^(x) — )■ 6{x). 
With the boundary conditions g{—L) = g{L) and differentiability on S, the one-particle 
Hamiltonian is clearly Hermitian under the inner product (■, ■) (I2.15p . Hence, by similar 
arguments as those of Appendix O, it gives rise, after completion, to a Hermitian operator 
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on the full Hilbert space L^(§) © C. The eigenvector equations fl2.10p for the even sector, 
with eigenvalue p, are replaced by: 



id^gpix) + Thp(f^''\x) = pgp{x) 



(2.26) 



using {gi,g2) ■= {gi © 0,5(2 © 0). Let ^p(x) := e *^^5'p(x) and (p^p\x) := e ^^^(/^['''(x). Eqs. 
f l2.26p then become idx^p{x) = Thip^p\x) and r {'^p'\C,p) = (p — e) hp, and their solutions 



are: 



gp{x) = -%T hp [$W (x) + Cp] e'P" with ${ 



where the integration constant is 



-p ■- 



(2.27) 



The inner product in Cp can re- written, using integration by part, as 



1 



dx^'';\xfe^'^\ 



(2.28) 



la 



The periodicity relation gp{—L) = gp{L) imposes the quantification condition: 

Note that |1 + <^^p\l+)/cp\ = 1 because 23fJe (</?jr^, 4"') = l4"'(4^)P- The quantity 1 + 
$|r'(/^)/cp is the phase shift incurred through the spread impurity. 

The large-L limit can then easily be taken, whereby the distribution of p values is 
uniform. With appropriate normalisation, we find, for the full operator solution. 



gp{x) 



1 + 



ipx h ^ _ 
TCp 



(2.29) 



and 



1pe{x,t) 

i)o{x,t) 

—IT d{t) 



— gp{x)e 'P^Qp, 



27r 
dp 



Jp{x-t) ^ 



2tt 
dp 



c„^e-'^Up. 



(2.30) 



The canonical anti-commutation relations give again f l2.24p . Note that the orthonormality 
relation 

dx gp{x)gpf{x) + hphp/ = 2n6{p — p') (2-31) 
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follows from the eigenvalue problem and the chosen normalisation, using the fact that 
|1 + $|r^(/^)/cp| = 1. This in particular gives the modes in terms of the fermion fields: 

/"°° dx i 
Up = -^gp{x)iJe{x) —d 

J -oo V 27r TCp 

/"°° dr 

bp = / e-'^^Mx). (2.32) 



The inverse orthonormality relations are simply obtained from fl2.3ip and from the fact 
that Qp © hp form a complete basis on L^(]R) © C: 

J dp gp{x)gp{x') = 2ti5{x — x'), j dphphp = l, J dpgp(x)hp = 0. (2.33) 

The Hamiltonian takes the same form as before, as well as the Hilbert space, in terms of 
the new operators and bp. 

The discontinuity condition (12.221) does not completely disappear: it is simply smoothed 
out. It can be written in the form 

MO - Mia) J ^ {e'''" - e'^'^) apC-^P' - zr j ^ dyv^'^\y)d{y - /+), (2.34) 

valid in any finite-energy state. Hence, there is again hybridisation. 

The limit a — ?■ 0"'' can be taken, independently of the limit L — )■ oo (both limits 
commute). In the limit a — 0"*", we have ^p^p\x) — > 6{x) and $p"'(x) ©(a^), the 
Heaviside function. Further, since \Ia\ and J+ — > 0, expression fl2.28p makes it 
clear that (y^jT', ^jf') — ?■ 1/2. In this limit of vanishing cut-off, the eigenfunctions Qp © 
/ip, with the normalisation chosen above, simplify to Cp © iwp/r, hence we recover the 
solution obtained from the resolved impurity above. Also, in this limit, the regularised 
discontinuity condition (12.341) goes to the discontinuity (12.221) . 

In summary, the regularised version (I2.25P of the Hamiltonian justifies the prescription 
consisting in defining g(fi) in Eq. (I2.10p as the symmetric sum (^'(O"'") -|- g{Q~))/2. The 
advantage of the regularised version (I2.25P is that it allows us to unambiguously define 
all averages. 



2.2.3 The case r = 

In the case r = (the "free" case), the solution to the equations of motion can be 
obtained similarly, with the principal difference that there is no hybridisation. Hence, 
for the fermion operators ijJe,ipo, we simply have to specialise (I2.17P or (I2.23P to r = 0, 
whereby both ipe and ipo take the same form; for the impurity operator, we simply have 
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d{t) = e'^^d. That is, 



^P'rix,t) = I ^e^f(-*)6p, (2.35) 

with anti-commutation relations fl2.24p for dp and bp, as well as {d\d} = 1, other anti- 
commutators vanishing. The lack of hybridisation is the fact that the corresponding one- 
particle solution uses a basis of functions in L^(M) © C where the subspaces and C 
are explicitly separated. Recall that at t = 0, flOBj) (or IHAT} ). flOOj) and f l235|) are just 
different parametrisations of the same fundamental operators: ip^Jl'^{x, 0) = ■ipe,o{x, 0) = 
ipe,o{x) and d^'^%0) = d{0) = d. 



2.3 The large-time limit and the stationary- state density matrix 
2.3.1 Description of the density matrices 

The initial density matrix pa ® pth of (12. 7p is most naturally constructed using mode 
operators of the r = solution: the Hamiltonians Hj and the number operators A'^- are 
explicitly diagonalised. More precisely, it takes the form 



Pth = exp 



dp [ f3p{alap + hlhp) H ^ (a^fep + h^dp^ 



/{Z,Z2). (2.36) 



Although the meaning of this density matrix as an operator on a Hilbert space may be 
subtle to clarify, we will adopt a simpler viewpoint, understanding density matrices as 
giving rise to linear functionals on an appropriate space of operators. From this viewpoint, 
the density matrix is an object - a measure - that tells us how to associate averages to 
operators. We do so by using averages of products of mode operators, which are calcu- 
lated formally using the cyclic property of the trace and the canonical anti-commutation 
relations. Hence, we define 

Tr(pth alpaq) = - 
Tr(pthSj6g) = - 

TT{pthblbg) = 2 (^1 + g(/3+vy/2)p-y/2 + ^ _^ ^(i3-w/2)p+v/2 j (2.37) 

along with application of Wick's theorem for averages of products of more operators. 

The non-equilibrium steady-state density matrix pstat given in (12. 8 p is defined similarly; 
it is expressed most easily using the r 7^ mode operators Op and bp (both in the resolved 



1 - 


^ g(/3+Vy/2)p-y/2 


H 


_ g(/3-vy/2)p^ 


-V/2 




1 




1 




1 - 


^ ^(l3+W/2)p-V/2 


n 


^ g(/3-iy/2)pH 


-V/2 




1 




1 
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regularisation and in the spreading regularisation). Indeed, these operators diagonahse 
the interacting Hamiltonian if, and a solution to the equations fl2.9p defining the operators 
Y and X is 

y = \j dp {a% + 6jap), X = ^jdpp (ajfep + fojop). (2.38) 

Again, using canonical anti-commutation relations and cyclic property, the functional 
evaluating averages under the density matrix pstat can be defined. Of course, by construc- 
tion we find the same formulae as in (12.371) : 

Tr(pstat alag) = Tr(pth aj^s) 
Tr(pstat aj&<?) = Tr(pthaj6g) 

TT{pst.tb%) = TriptJlk) (2.39) 
along with application of Wick's theorem for products of more operators. 

2.3.2 Uniqueness of the (generalised) Hershfield operators 

In order to show uniqueness of the steady-state density matrix, constructed using the 
operators X and Y as defined around (12.91) . we only have to show that for any operator 
W such that [IV, if] = and [IV, ?/;e(a:)] = [IV, '?/'o(x)] = for all x small enough, we 
have W (X 1. It is sufficient to show that [IV, a^] = [IV, 6p] = for all p G M, where 
ttp and bp are the modes of the solution (12.231) or (I2.30p for the fermion operators in the 
resolved or spread impurity regularisation. Thanks to [IV, H] = 0, we have [IV, ipd^, t)] = 
[W,ipo{x,t)] = for X small enough, and for all t G M. Hence, in both regularisations, we 
find (with y = X — t) 

Fourier transforming in y gives the result. This indeed shows that the definition given 
around (12. 9p is a good definition. Note that we only need to have one value of x, small 
enough, where (12.91) holds. Note also that this proof immediately holds as well in inter- 
acting impurity models (with free-fermion baths), because for x small enough, the modes 
of the fermion operators are the asymptotic state creation operators [5]. 

2.3.3 Existence and description of the large-time, large-L limit 

We now show that the limit (12. 7p exists in a weak sense. 

Proposition 2.1 Define a finitely supported operator G as a finite linear combination of 
finite products of the form A HLi ^i(^j); withOi G {■ipe,'ipo,ipl,4'l} and A G {l,d,d'^ ,d^d}. 
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Then the limit 1^2. T\) exists in a weak sense, and is described by the non- equilibrium steady- 



state density matrix Pstat- 

lim lim Tr (Uo.to (Pd ® PthWl^^c) = Tr (p^.^.G) . (2.40) 

to— oo L— >oo \ / 

We show these statements both in the resolved-impurity regular isat ion, and in the spread- 
impurity regularisation. A nontrivial aspect is that the discontinuity equation f l2.22p holds 
in the resolved-impurity case in the trace above at to < but not at to = 0; from the 
spread-impurity calculation, we also show that the discontinuity equation arises according 
to two time scales, one being the spreading length \Ia\, where the discontinuity is obtained 
up to terms of order \Ia\, the other being the time scale characterising the approach to 
the steady state, of the order of 2/r^ when \Ia\ is small. In particular, for the resolved 
impurity, the discontinuity appears immediately. 

• Resolved impurity: 

The proof relies on writing the solution to the equations of motion with r 7^ in terms 
of local operators. It is a simple matter to see that the equations 

-idtipo{x,t) = idr,4jo{x,t) 
-idtipe{x,t) = id^il)e{x,t) - T5{x)d{t) 

-idtd{t) = -Tipe{0,t) - ed{t), (2.41) 

where ipe{0,t) = {ipe{0~^,t) + ipe{0~,t)/2 for any t > 0, is solved as follows: 

1po{x,t) = 1po{x-t) 

Mx,t) = ipeix - t) - iTQ{x)Q{t - x) ( -tr [ dx' ipe{x' + X - t)e-'"'' + e-^^'-^'U 







t 

bx' I -bt 



d{t) = -IT dx'iPeix' -t)e-^-^ +e-^'d (2.42) 



for < t < 2L and where b := r^/2 -\-ie. Note that the discontinuity relation ( 12.22^ holds 
exactly for all t > 0. In the expression 

Tr (Uo;t,iPd ® Pth)f/oU g) = Tr {{pd ® Pth) G(-to)) (2.43) 

we have all local operators in G'(— to) evolved at time t = —to. This is evaluated using 
the solution fl2.42p . and using the form (12.361) of the initial thermal density matrix along 
with the expressions (12.351) . Note that the impurity operator part simply factorises, and 
that the electronic part can be evaluated using formulas (12.371) . The result is a linear 
combination of multiple integrals over positions x' of finite products of local operators 
ipe{'),'^o{'),d. The integrands of the resulting integrals may have poles where the positions 
of the local operators collide. The prescription for the evaluation of the integrals is to 
deform the integration contour in such a way that the complex positions of operators have 
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negative imaginary parts that are in increasing (towards the positive imaginary direction) 
order from the left to the right. This prescription guarantees the correct ordering of the 
product of operators and that the local-operator solution is in agreement with the mode 
solution discussed above. 

The limit L — )■ oo can easily be taken, and the only effect is to make the explicit 
local-operator solution f l2.42p valid for all of < t < oo. Then, the limit to ^ —oo can 
then also be taken. For any given G, there is a (negative) value of to large enough such 
that the factor Q[t — x) in the solution fl2.42p for ipeix, t) is, when the solution is put into 
G(— to), always 1. Also, since correlation functions of local operators in the density matrix 
Pd®Pt\i do not grow at large separations (in fact, go to constants), and since the integrand 
in the solution (12.42^ has an exponentially decaying factor e~^^ , the large-time limit of 
the multiple integrals exist. Since, moreover, the operator d in the solution f l2.42p always 
come with an exponentially decaying factor e~^*, all terms containing it on the right-hand 
side of fl2.42p can be omitted in the large-time limit. Hence, in pd ® Pth, we may omit 
Pel- As a result, since the ensuing correlation functions only contain the local operators 
ipe{x — t), ipe{x' + x — t) and ipe{x' — t) (at various positions x and with integrations under 
x'), by space-translation invariance of the density matrix Pth and of the contour-shift 
prescription, we may omit the term —t in the arguments of all local operators in (12.420 . 
Hence, the limit L — )■ oo and to ~^ —oo of the average of G(— to) can be obtained from 
the modified, infinite-time local-operator solution 

^^o(a;, oo) ipo{x) 

/■oo 

ipe{x,oo) ipeix) - t'^Q{x) / dx' %Ij{x' + x)e~^'^' 

Jo 

POO 

d{oo) ^ -ir / dx' iJeix')e-'"'' . (2.44) 







Replacing ipei'), V'o(-) by their r = mode expression (I2.35p . we find 

M^,oo) ^ [ -^e'^^bp, (2.45) 
J V 27r 



ilJe{x,oo) J -^=ep{x)a. 
—iTd{t) J 



dp 



We see that the right-hand sides are exactly the expressions f l2.23p for the non-zero-r 
solution. Since these must be evaluated against pth, and since the evaluation of r = 
modes against pth, (I2.37p . gives the same answer as that of non-zero-r modes against the 
non-equilibrium density matrix Pstat, (I2.39p . this proves that the limit L — )• oo, to — )■ — oo, 
in this order, of pth gives pstat when evaluated against finitely supported operators. 
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Spread impurity: 

The equations of motion now take tlie form 

- idtipo{x,t) 
-idti^e{x,t) 

—idtd{t) 



idxipo{,x,t) 

id^Mx,t)-Tip^''\x)d{t) 

dxip^°'\x)tlJeix,t) - edit). 



— T 



(2.46) 



la 



A solution is obtained through 

^Peix,t) = ij,{x-t)-iT [ dy ip^''\y)Q{x - y)Q{t - X + y)d{t - X + y) (2.47) 

J la 

for < t < 2L, where d{t) satisfy the integro-differential equation 



-idtd{t) 



-T 



dx' (^["^l (x')^Ae(x' -t) + iT^ / dx' (f)^"^ {x')d{t - x') - edit) 



with 



rfy^W(x' + 2/)^W(2/). 



(2.48) 



la 



This integro-differential equation can be solved recursively as a power series in r; the power 
series has an infinite radius of convergence, for any finite t and when evaluated inside 
averages with other local operators. The structure of this solution is not so important, 
except for the fact that the limit L — )■ oo can, as before, be taken easily, the only effect 
being that the solution becomes valid for all t > 0. 

Recall that from Eq. (12.431) we first have to evaluate expectations specified by the 
initial density matrix Pd® Pth of local operators G(— to) evolved at time t = —to and, then 
to take the limit to ~^ —oo. So we have to look at the behavior of local operators at large 
time t — )• +00. 

It is of course a simple matter to obtain a more explicit form of the L = oo solution, 
using the time-evolved operators ipe{x,t) and d{t) in (I2.30p along with the expressions 
(I2.32P of the mode operators in terms of the initial conditions: 



'^e{x,t) = — 

Ztx 

d{t) 



2txt 



dx' 
dx' 



dpgp{x')gp{x)e 



i)e{x') 



2txt 



d 



dp ■ 



\x' 



-ipt 



^e(x') 



1 

2^ 



dp 



-ipt 



d. (2.49) 



The explicit solutions f l2.29p and fl2.27p make it clear that as lm{p) — t- — oo, the expressions 
gp{x) / Cpe~^^^ and l/|cppe~*P* vanish exponentially if t > |/^| + + Hence, in the 
solutions above, the integrals in the part proportional to the operator d can be evaluated 
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by such a contour deformation. The result IS, cLS cL function of time t, exponentially 
decaying. Hence, at large time, the part in d vanishes, so that we may restrict ourselves 
to the part containg ipei^') (and in pth? we may omit p^). 

Note that the leading exponential of the large-time decay is of the form e"^"* where 
the (positive) real part of ba is the (absolute value of the) imaginary part of the position 
of the first zero of Cp away from the real axis. The expression f l2.27p indicates that, for a 
small enough, this will be near to r^/2 (more precisely, lima^o^a = t^/2 + ie). 

In order to correctly evaluate the part containing ipe{x'), we follow the technique used 
in the resolved-impurity case, and shift x' to x'+t, without changing the result of averages. 
Hence, we have 



d{t) ^ 



2txt 



dx' 
dx' 



dpgpix')gp{x)e 



-ipt 



Ipeix' + t) 



as t — oo. The limit of large t can now be evaluated by replacing ipe{x) by its free- mode 
expansion, 

1 



d{t) ^ 



(27r)3/2 



dq 
dq 



dp I dx' gp+g{x')g. 



p+q 



X e 



iqx' —ipt 



(27r)3/2r 

We evaluate the integrals over x' as follows: 



dp I dx 



9p+q{X ) ^iqx'-ipt 
Cp+q 



(2.50) 



dx' gp+q{x')e 



/\Jqx 




dx' -gp+q{x')e"^^' 

p-ipla 



-ip + 0+ 



+ 1 + 



$1 (/ + 

p+q\ a 
Cp+q 



-ipla 



ip + 0^ 



By the same argument as above, since gp oc 1/Cp, the integration over p of the first integral, 
occuring when it is put back into the expressions in the square brackets in fl2.50p . gives a 
quantity that vanishes exponentially as t — oo (again, as e"^"*), hence can be neglected. 
Further, the last term is subject to the same argument, since the additional explicit pole 
is at p = iO~^, hence is not taken as lm{p) — > — oo. The only non- vanishing term as t — ?■ oo 
is the one occurring from the pole a.t p = —iO~^ in the second term. The result is 

dp 



V'e(a;,oo) h-^ 
—iTd{oo) I— i- 



s/2Tr 
dp 

dp 
v/2^ 



Jp{x-t) ^ 



9p[x)e 



-ipt ~ 



"p^ 



(2.51) 
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By the same argument as in the previous paragraph, this reproduces (12.301) if the density 
matrix pth is replaced by Pstat, hence shows the proposition. 

Finally, we see from the solution f l2.47p that the discontinuity condition f l2.34p is 
replaced by 

^e(/+,t)-^e(/-,^)=^e(/+-t)-^e(/--t)-^r / dy ^^'^\y)Q{t - 1^ + y)d{t - 1^ + v) . 

J la 

Hence, at t = we recover ipeila) ~ V'e(-^a") right-hand side, whereas at t > \Ia\ 

we recover something that looks similar to fl2.34p . except for the first two terms. The 
discussion above, showing that as t — )■ oo we can re-interpret the free modes dp as the 
modes a^, implies that the first two terms tend to the first integral of fl2.34p . Hence, we 
see two time scales: the time scale \Ia\, after which the main part of the spread-impurity 
discontinuity is reached up to an error of order \Ia\, and the time scale 1/Re(fea) ~ 2/r^ 
necessary for this error to go to zero. 



3 Transferred charges and the Levitov-Lesovik for- 
mula 

We are going to present the proof of the Levitov-Lesovik formula for the p.d.f. of the 
transferred charges in the case of the resolved impurity. It is clear that an analogous 
proof can be done in the case of a spread impurity: the logic will be the same but the 
notations and the computations will be much more cumbersome (without much gain). 



3.1 Definitions and general comments 

As explained in the introduction, we aim at describing the charge A(t) = Q{t) — Q 
transferred during a large time interval t in the stationary regime. The charge is chosen 
to be that in the first reservoir, so that 

Q = Ni = J dx iljl{x)ipi{x) 

and, as usual, Q(t) = Ul^QUtfl- 

Let the system be prepared at time to, and denote the density matrix at time zero by 
Po := Uo;t^{pd <S) PthjUltg- The first definition flLip of the charge p.d.f. gives 

oo 

Ptoi^,t)= e^^" Tr (p,,+,f/iP,,poP,„f//no+.) • 

q,qo = ~co 

The second definition of the measured transferred charge would be similar but with only 
one projector Pg^. See Eq. ( 11. 2p . 
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The first sum can done using the relation J2q f{(l)Pq = f{Q)- Since Q has an integer 

spectrum, the second sum can be dealt [9] using the formula J^^ dfie'^^^^~'^^ = 27rPg, so 
that 



p^^(A,t) = / " ^ Tr (p„e-*("/2+;.)QgaQ(t)g-i(A/2-;.)QN _ ^3 52) 

Jo 27r 

The second definition for the generating function yields a similar expression for Ptg{X,t), 
the only difference being that /i is set to inside the trace (so that the //-integration does 
not matter). 

In order to reach the stationary regime, we need to take the limit L ^ oo then 
to — )■ — oo of Ptf,{X,t). Clearly, neither Q nor Q(t) for any t is finitely supported; hence, 
their averages do not admit a large L limit. Physically, this is because the reservoir carries 
an infinite charge at L — )■ oo. But the operator 0A,/i(t) defined by 

gS^.^W _ g-i(A/2+M)QgiAQ(t)g-i(A/2-M)Q (3.53) 

is finitely supportecH. 

Indeed, first, the charge difference A(t) := Q(t) — Q is finitely supported, with spatial 
extension of size vpt, because the term in Q{t) diverging for L large is time independent. 
More precisely, we find 



A(t) = [ dt'i[H,Q]{t') = [ dt'Ji{t'). (3.54) 
Jo Jo 

where 

Ji(t) = Ty/2lm [^1(0, t^d] (3.55) 

and we see that A(t) is supported on the segment [—VFt,0]. The large-L limit of this 
operator can be taken from this expression. In terms of modes, the operator A{t) is 

(the restriction on p 7^ /c in the sum is automatically ensured by = 1) with 

ipt _ 1 

MP) = (3.57) 

One can also verify directly from this expression that the large-L limit exists. However, 
the calculation is more subtle. In particular, the naive regularisation of the large L limit 



■^To be more precise, it lies in a completion of the space of finitely-supported operators, in the sense 
that it involves a converging infinite series of operators uniformally supported on a fixed finite region. 
However, although every finite sum converges towards the steady state according to Proposition [^Hl we 
will not go into the details of the convergence of this infinite series towards the steady state, as this is 
beyond the scope of the present paper. 
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of A(t) which would consists in transforming the sum into an integral and regularising 
the pole in l/{p — k) coming from 6t{p — k) with an +iO~^ prescription, l/{p — k) i— )■ 
l/{p — k ± iO~^), is incorrect. 

Second, we write Q{t) = Q + A{t) and use the Baker- Campbell-Hausdorff formula to 
write the product of exponentials as an exponential of multiple commutators, 

0A,/x(^) = iXA{t) + commutators between Q and A(t). 

The commutators are all operators supported on the segment [—vpt, 0] because A(t) is, 
and because Q is an integration of a local density. Hence, by Proposition 12. the average 
of any function of Q\^fj,(t) admits a stationary limit, and 

Pstat(A,t) := lim lim PtoiX,t) = / ^ Tr (p.tate®^-^) . (3.58) 

to— oo L—^oo J ZTX 

Finally, the universal behaviour, where the Levitov-Lesovik formula holds, emerges in 
the large t limit. This is the limit tr^ ^ vp since vpjr^ is the damping length in the 
stationary regime. We want to evaluate the large deviation generating function F{X) for 
the charge statistics, defined by ( 11. 3p . 

Note that we shall deal only with the first definition (11. ip for the generating function, 
so that Eq. (I3.58P applies. But we shall later show that Tr {^pstat e®^'^ ^*-*) is actually /i- 
independent at large time, so that the large deviation function is identical for the two 
definitions. 



3.2 Expression for Oa,^(^) at large L and large i 

In order take the limit L — j- oo in the definition (I3.53p . we must perform the multiplications 
of the exponentials in terms of explicitly finitely supported operators. Since Q[t) is 
a bilinear expression in the fermion operators preserving fermion numbers, it lies in a 
representation space of the Lie algebra associated to the group of linear endomorphisms 
of the one-particle Hilbert space. Hence, 0A,^(i) also lies in this representation space. 
Let us denote by g, g(t) and Q\p.{fy the one-particle operators associated to Q, Q{t) and 
0A,/i(t), respectively. Clearly. 

Further, let 6t, S'^ and 5^ be the one-particle operators associated to A(t) = Q(t) — Q, to 
[Q,A(t)] and to [Q, [Q,A(t)]], respectively. These are all finitely supported operators. 

Proposition 3.1 The one-particle operator 9 \^^{t) corresponding to Q\^^{t) is given by 
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where 

cs = i sin A, C52 = —2 sin^ — , = 2 sin /x sin — , 

Proof. Notice that g is a projector, = q, because on the one-particle subspace, Q is 
the operator counting the number of electrons, hence on the one-particle Hilbert space, 
there is a basis that diagonalises with eigenvalues and 1. Similarly, q{t) = q + 6t is also 
a projector, q{t)^ = q(t). Hence, 

e^A.MW _ 1 = + y^q§^ + y^S^q + yiy2q5tq)y3 

where yi = e"*'^'^/^"'"'^^ — 1, y2 = e~*^'^/^~^^ — 1 and j/3 = e*'^ — 1. Using = q, we have 
further 5^' = 6t — 2q6tq + Stq. Using also (g + 6tY = q + 6t, we have q6t + Stq = St{l — St), 
so that 

2q6t = 6til-6t) + 6't 
2 6tq = 6t{l-6t)-6[ 
2q6tq = 6til-5t)-5':. 
Inserting this set of relations in the previous formula for e^^'^'^^\ we obtain Eq. f l3.59p . □ 



Expression (13.591) for the operator e^^-'"^*-' is explicitly in terms of local operators, hence 
the limit L — )■ 00 can be taken. More precisely, the matrix elements of 5t, 5[ and 5'1 may 
be calculated using expressions of A(t), [(5,A(t)] and [Q, [Q, A(t)]] as bilinears in the 
modes Op and hp of the large-L theory obtained from expressions (I3.54p and (I3.55P for 
A(t). Using (I2.42p . we may define two operators A{t) and B{t) and express them in the 
large L limit in terms of the fermionic fields: 

A(t) := lim [Q, ci(t)] = ^/ cixV^i(x-t)e-("'/2+M^ (3.6O) 

V2 Jo 

St(t):= lim[Q, 7^1(0, t)] = + i ^[Q . d{t)]\\ (3.61) 

L-i^oo L^oo \2v2 / 

They can also be expressed in terms of the modes and hp-. 



Bit) = ^ / ^ {i^P + 4)e-^^* - ^pe-(-^/2+-'^)*) (a, + bp) 

From dSSni) and dMH), we have [g, [Q,d{t)]] = -[Q,d{t)] and [Q, [Q, 7/;|(0, t)]] = [Q, 7/;|(0, t)], 
so that we find 

lim [Q, ijliO, t)d{t)] = B\t)d{t) + V^|(0, t)A{t) 
lim [g, [g, ^1(0, t)d{t)]] = B\t)d{t) + 2B\t)A{t) - ^1(0, t)A{t) (3.62) 

L— >oo 
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From this and the equivalent mode decompositions of d(t) and ipi{0,t), we obtain the 
mode decomposition of A(t) and of its commutators with Q. 

We need these expressions for t — > oo. In order to obtain the leading large-t limit, we 
can keep in the expressions for A{t) and B{t) only the purely oscillatory exponential, the 
other terms lead to exponentially decreasing factors. Hence, upon integration over t and 
after some algebra, we find the one-particle operators: 



with St{p) 



IT 



Q,-^ / dt' iljl{Q,t')d{t' 
v2 Jo 



Similarly, we have 



2wp Wp — Wp 



Wp - Wp 



-2Wr. 



0{e 



-^r2/2^ 



Q, 



Q,-^iyt'4io,t')dit') \p') ^ ^sip-p' 



Wp — Wp —2wp 
2wp Wp — Wp 



Adding the hermitian conjugates and putting everything together, we obtained the desired 
expressions: 



{p\S,\p'} = 


Ihp 






(pW) = 


Is.iP 


mv) = 





Wp 
Wp 



Wp - Wp 
















Wp - Wp 



Wp-Wp ^ , nf.-t'T-'n 



Wp - Wp 



+ 0(e- 



(3.63) 
(3.64) 
(3.65) 
(3.66) 



We can now put this into (13.530 to obtain: 

Proposition 3.2 In the large L, and then large t limit: 

A 



with 



M{p) = e'^"^ 



ia^i sin — 2 sin — 1 sin — h ia^. cos — 
y -rp 2 \ 2 2 



A 



(3.67) 



(3.68) 



where (pp is defined in /12.18\) . 



Proof. Using f l3.59p and the previous expressions f l3.63tl3l66|) . we write e^^'"^*-' as in 
Eq. (13:671) with 



M{p) :-- 



s+C-Wp + c+S-Wp ic+C-Wp — is+S-Wp 
ic+C-Wp — is+S-Wp c+S-Wp + s+C-Wp 
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where s± := sin ± ^) and c± := cos (| ± Some algebra also gives this in the form 

M(p) = i sin (j)p {a^ sin fi + Cy cos /i) — 2 sin^ ^ ^1 sin — + icx cos — ^ , 
which is equivalent to Eq. (13. 68 p . □ 



3.3 Determinant formula 

We want to evaluate the large deviation generating function F{X) for the charge statistics: 

F(A) = limrMogPstat(A,t). (3.69) 



i— >oo 



using fl3.58l) . The formula obtained above for the one-particle operator Ox^^{t) correpond- 
ing Q\^^{t) can be put to good use from standard free fermion technology. We rewrite 
determinant as follows. 

Proposition 3.3 We have 

det(l + ristat) 

where 9\^(t) and Ustat are the one-particle operators associated to Qx,^(t) and pstat respec- 
tively, and Gstat = i^n'^l t rnatrix of fermion two-point functions. 

Proof. Notice that Pstat only involves products of operators that are elements of the group 
of linear endomorphisms of the one-particle Hilbert space. This follows from the discussion 
above for e®^ ''*^*'', and from the fact that the operators H, Y and X are exponentiated 
in Pstat, Eq. (12. Sp . are bilinears in fermion operators that preserve the particule number. 
Further, recall the formula Tr(e^'^) = det(l + e™) where M is such a bilinear in fermions, 
and m is the associated one-particle operator. Here, Tr is the trace on the full Hilbert 
space, and det is the determinant on the one-particle Hilbert space. Clearly, by the 
group representation property, this formula extends to product of similar exponentials 
of bilinear, e.g. Tr(e^'^e^) = det(l + e"^e") (with obvious notation). Hence, we obtain 
formula f l3T0|l . □ 



3.4 Levitov-Lesovik Formula 



Since for t — oo, the one-particle matrix Ox^^{t) becomes block- diagonal in the momentum 
basis as in (I3.67p . with a leading term proportional to 5t{j) — p'), then as t — oo we find: 



log Tr (e®^-Wpstat) =t / ^ log 



det(l + n{p){e^^'^^f^ - 1))] + 0(1) 



(3.71) 
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where the 2x2 matrices Ox^^{p) and n{p) are the blocks of momentum p of the one-particle 
matrices associated to the large time asymptotic of 0A,^(i) and to Gstat? respectively. The 
proof that the term St{p—p') gives rise to the linear dependence at large t upon evaluating 
the trace (or the determinant) is presented in Appendix |A] (Equation OA. 741) ). 

The matrices involved in fl3.7ip are, explicitely, 

ge^,M(p) _ 1 = sin ^ M(p), n{p) = ^ [(ni(p) + n2(p))l + {ni{p) - n2ip))a,^] 

with nj{p) the thermal occupation numbers of the reservoirs, defined in (11. 5p . Since h{p) 
commutes with a^, the form of M{p) in f l3.68p makes it clear that the leading term in 
f l3.7ip is independent of Hence, the two different ways (II. ip and (II. 2p of defining the 
measure of transferred charges give identical large t limit, and the large-deviation function 
is thus given by: 

FiX) = /^log [det(l + sin^M(p)n(p))]. (3.72) 
This shows the following: 

Proposition 3.4 The large- deviation function F{X), both as defined in ( li.ij) and in ( li.i^j) . 

is given by the Levitok-Lesovik formula [T^. 



Proof. This follows from (I3.72p . taking into account the dispersion relation uj = p and the 
transmission coefficient (I2.19p . □ 



4 Conclusion 



We have provided a description of the resonant-level model and its non-equilibrium steady 
state based on Hershfield's non-equilibrium density matrix, showing that this density 
matrix is obtained for averages of local (more precisely, finitely-supported) observables, 
and clarifying aspects related to the point-like impurity. We have then given a precise 
derivation of the Levitov-Lesovik formula based on local obsevables, making it clear how 
it enters into the framework of Hershfield's non-equilibrium density matrix. A direct 
extension of this work is an extension the derivation to a perturbative analysis of the 
(deformed) Lesovik-Levitov formula for the interacting resonant-level model. This is 
natural because of the local-operator nature of the derivation, well suited to perturbative 
analysis. On the other hand, it would be very interesting to study the full-counting 
statistics in interacting integrable models - for this, the scattering approach, taken by 
many authors for deriving the Levitov-Lesovik formula, would probably be more useful. 
Hence, an important step is to understand the relationship between the correct local- 
observable description of the steady-state and the scattering approach, in order to have 
a precise scattering-state derivation. We hope to come back to these questions in future 
works. 
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A Large-t behaviour of the regularised delta-function 



We wish to evaluate the right-hand side of the last equation of (13.701) in the large-L, then 
large-t limit. According to Proposition 13. 2[ this has the form 

det(l + AB{t)) 

where the determinant is evaluated on the one-particle Hilbert space. This Hilbert space is 
a tensor product of a space with a basis enumerated by a continuum momentum variable, 
and a finite-dimensional ( "internal" ) space. The operators have the form 

{p\A\p') = 6{p-p')A^, {p\B{t)\p') = 6t{p-p')B, + 0{e-'^'/') 

where Ap and Bp are matrices acting on the internal space. 

The determinant on a continuous-basis Hilbert space can be defined via the formula 

logdet(l + AB{t)) = Trlog(l + AB{t)) = V ^-^ Tr [{AB{t))''] 

k=i ^ 

where the trace is defined via an integration over the basis states, 

Tr [{AB{t))'] = tr dp {p\{AB{t))'\p)^ 

(with tr the trace over the internal space). Hence, we have to evaluate 

k 

I dp{p\{AB{t)f\p)= I dp^---dpu\[Ap^Bp^6t{p,-pi+^) 

(with Pk+i = Pi and the product being ordered) in the limit where t — )■ oo. The most 
direct way is to use Fourier transforms for the momentum dependence (with coefficients 
being matrices on the internal space). We will show that, in the sense of distributions 



over Q!jS, 



k 

j dp,--- dp, e'^U ".f. \[ 6tip, - p,+i) =t [y + ^(1) (^-73) 

.7 = 1 ^ 



as t — oo. This implies that 



logdet(l + AB{t)) =tj^ det(l + ApBp) + 0(1) (A.74) 
where the determinant on the right-hand side is on the internal space. 
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In order to show (lA.73p . we use fl3.57p . in the form 

Hence, the left-hand side of (I A. 731) becomes 

^ 1^ dt[ ■ ■ ■ dti j dp,... dp, e^Ej^i^'.fe-*.-+-.). 

Evaluating the integrals over pjs in terms of delta-functions, we find that all tj-integrals 
are cancelled except for one: 




where A = max(Q;2, 02 + «3, • • • , J2j=2 '^i) ^ ~ niin(a2, ^2 + «3, • • • , J2j=2 '^i)- "^^^ 
large-t behaviour gives the right-hand side of (]A.73p . 



B Current fluctuations 

This Appendix is dedicated to show how the statistics of transferred charges is very 
different from that of the integrated current. As we shall see the latter is non-universal, 
depending on a high energy cut-off. Although the mean of the integrated current is finite 
(because it is equal to the mean of the transferred charge), its p.d.f. is made of three 
peaks, one centred at zero current, and two others centred at opposite positions diverging 
with the high energy cut-off. Thus measures of the local current on the impurity are 
highly fluctuating variables. 

We shall only deal with the resolved impurity case. The integrated current over a time 
interval t is defined by Xt = ds Ig with 1^ = r'^m ['^1(0, s)d{s)~\ . We are interested in its 
p.d.f. so we look for the generating function of its cumulants defined as the expectation 
of e*""^' in the stationary state: 

F™t(a) := Tr(pstate^^^') = Tr(pstat e^*^^°) , 

where we used the stationary property in the last equation. Since Iq is bilinear in the 
fermion operators, the trace can be evaluated using a formula analogous to Eq. f l3.70p : 

Fcurrent(a) = dct [l + Gstat(e**"'" - 1)] , 

withTo the one-particle operator associated to Iq and Gstat the matrix of fermion two-point 
functions. 

We skip the details of the computation of this determinant which needs to be reg- 
ularised. Let A be a large momentum cut-off: \p\ < A. The result has the following 
form: 

i^™t(a) =Ao + e+^*"v/^^ + e'^'^V^^ A., (B.75) 
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with amplitudes Aq and A±. given by: 

A+ + A^ = J [(V^IVo)a - 2{ijliJo)A {Sd)A + 2(V'Ic/)A(dVo)A] + {d^d)A 
Aq = 1-A+-A. 

where (■ ■ ■ )a are the regularised stationary expectations. So the peaks are centred at 

C The Hilbert space L'^{S) C 

Consider the functions ep{x) defined in fl2.12p . These functions have the properties 



ep{x + L) = ep{x - L) , ep{ 



~Xj CpyX jVp^ 



as well as 

^ ( 2L {p = p') 

dx epi{x)ep{x) = < WpiWp (C.76) 



{P 7^ P') 



where Vp and Wp are define in fl2.13p . 

Consider the linear space S of sequences (a^ : e^*^^ = Vp) (which we will also denote 
simply by (cp : p)) satisfying the asymptotic condition ap = o{p~^) as \p\ — )■ oo. Consider 
also the subspace J-" of L'^{S) © C, with elements which we will denote g(Bh, characterised 
by the fact that the wave functions g are continuous everywhere except for a finite jump 
at 0, along with the linear relation 5'(0~'") — g{0^) = —irh (here we do not need to specify 
the value of at 0). We now show that 

9{^) = ^^Cp(3^)'2p! h = — ^^WpOp (C.77) 
p p 

is a one-to-one map M : 5 — )■ J-", (op : p) \^ g ® h. 

To show that M maps S into J-", we evaluate, for 77 > small and x 7^ 0, 



I (x -|- 77) — (^(x — ?7)| < 2 I ap I I sin p?7 1 < 2//"'"/^ |ap| -|- 2 



Op I 



The right-hand side tends to as — )■ 0, since the sum of \ap\ is convergent. Further, the 
case X = gives 

9{ti) - 9{-v) = J^^Vpe'P'^ - e-'P'')ap = J^WpapC^^'' + 2i J] Op sinpr/. 

p p p 

28 



By the same arguments as before, the last sum on the right-hand side of the last equation 
tends to 0, and by convergence of the sum defining h in flC.77p . the first sum tends to 
—irh. To show that the map is one-to-one, we simply inverse the equations (IC.77p . By 
interchanging sum and integral (which can be done by uniform convergence) and using 
( 1C.76I) . we see that 



^ j dxg{x)ep{x) = ap + 



1 Wpiirh + Wpap)) 



2L T 



2 



hence ^ 

«P = — ( / dxg{x)ep{x) ^ ) . (C.78) 

2L + ^— 1^- \J-L 



Consider the inner product on J-" induced by the standard one on L^(]R) © C: 

{g®h,g' ®h') = j dx g{x)g'{x) + hh' . (C.79) 

This in turn induces an inner product on S (which we will likewise denote by (■,■)), 
obtained by the mapping flC.77p via ((a^ : p), (a^ : p)) = (M(ap : p),M(ap '■ p))- Using 
( 1076]) . we find 

{{ap : p), {a'p : p)) = ^ ^21 + a^- (C.80) 



We see that (1C.78I) can be written as 

1 



a. 



p 

2L 



-{g ®h,ep®iWp/T). {CM) 



Likewise, by inserting this into f lC.77p . we find that the vectors and u (the former 
having a wave-function part that is in fact a distribution) defined by 

^ ep{x){ep®iWp/T) ^ Wp{ep®Wp/T) 

^-'■=1^ or , K|2 ' ^-=2^ or , l»p|2 (^'^2) 



\wp\2 

p 

have the properties 



2L + ^ Y 



g{x) = {g®h,6^), h={g®h,uj). (C.83) 

Since the map M is one-to-one, we see, by the definitions of the inner products, that it 
is an isomorphism of the inner product spaces S and M(iS) C The completions of two 
isomorphic inner product spaces are isomorphic Hilbert spaces, and the map M extends to 
an isomorphism of these completions. The completion of S is the Hilbert space of square- 
summable sequences. We now show that the completion of M(iS) is the Hilbert space 
L^(§) © C. The main argument is to show that this completion must contain a direct 
sum of functions with an arbitrary discontinuity at 0, plus C; that is, that the linear 
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relation between g and h disappears. The completion of this direct-sum space under the 
direct-sum norm flC.79P is L'^{§>) ©C. To show the former assertion, we may form Cauchy 



sequences of functions in J-" such that the limit is a function where the discontinuity at 
is changed from h to another value - such sequences are straightforward to find. However, 
since we haven't shown that M(5) = J-", it is not evident that the sequence will lie in M(iS) 
without explicitly exhibiting an example. Consider the mapping ( ]C.77I) . with sequences 



(op : p) that are of the form Qp = c/p + o{p) as p ^ ±00, where c is some complex 
numbers. These sequences are in the completion of S, and we may show that they give 
rise to functions g whose discontinuity at can be any complex number. Indeed, in this 
case, the discontinuity at can be calculated as follows: 



p p ^ 



/sin p 
dp = —irh + 27rzco (C.84) 
p 

as ?7 — )■ 0"''. Likewise, a similar calculation shows that we may put any discontinuities 
at any finite number of points x. The completion of a space of functions continuous 
everywhere except for an arbitrary number of arbitrary discontinuities is L^(§). This 
shows that the completion of M(iS) is Iv^(S) © C. 

Finally, it is simple to see that every one-term sequence maps under ( 1C.77P to a solution 
to the equations of motion fl2.10p with finite energy E = p. Hence, the completion of the 
space of finite-energy solutions is indeed the Hilbert space i^^(S) © C (with the standard 
norm) . 
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